Dislocation-induced superfluidity in a model supersolid 
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Motivated by recent experiments on the supersolid behavior of *He, we study the effect of an 
edge dislocation in promoting superfluidity in a Bose crystal. Using Landau theory, we couple the 
elastic strain field of the dislocation to the superfluid density, and use a linear analysis to show that 
superfluidity nucleates on the dislocation before occurring in the bulk of the solid. Moving beyond 
the linear analysis, we develop a systematic perturbation theory in the weakly nonlinear regime, and 
use this method to integrate out transverse degrees of freedom and derive a one-dimensional Landau 
equation for the superfluid order parameter. We then extend our analysis to a network of dislocation 
lines, and derive an XY model for the dislocation network by integrating over fluctuations in the 
order parameter. Our results show that the ordering temperature for the network has a sensitive 
dependence on the dislocation density, consistent with numerous experiments that find a clear 
connection between the sample quality and the supersolid response. 

PACS numbers: 67.80.B- , 67.80.bd, 61.72.Lk 
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I. INTRODUCTION 

In 2004 Kim and Chan^'^ discovered a remarkable 
anomaly in the response of ^He crystals to ac rotation. 
Their results, together with subsequent experiments, '^^^^ 
have a compelling interpretation as the nonclassical rota- 
tional inertia (NCRI) expected for the elusive supersolid 
phase of matter. ^'^^^^ However, there are puzzles~for in- 
stance, the sensitive dependence of the Kim-Chan effect 
on sample preparation and quality.'^ Indeed, theoretical 
studies suggest that a pristine '^He crystal is "insulating," 
in that it does not support off-diagonal long range order 
(Bose condensation).^^ The emerging consensus is that 
sample defects, likely in the form of dislocations, play 
an important role in explaining disparate experimental 
results. Most existing theoretical works have focused 
on single dislocations; for instance, quantum Monte Carlo 
studies"'^^ have shown the existence of superfluidity along 
the core of screw dislocations in a model ^He crystal. In 
a phenomcnological approach. Landau models show that 
superfluidity nucleates first on edge dislocations in a su- 
perconductor or a Bose solid, ^^"^^ and Shevchenko^^'^^, 
Toner^^, Bouchaud and Biroli^^ have worked out some of 
the properties of a network of dislocations in a quantum 
crystal. In this paper we continue and extend the work 
of Shevchenko and of Toner by systematically deriving 
an effective random-bond XY model for a network of 
superfluid dislocations. 

This paper is organized as follows. We first introduce 
our model for superfluidity in the presence of a single edge 
dislocation, by coupling the elastic strain field due to the 
dislocation to the superfluid order parameter. Then, 
using a linear stability analysis, we show that superflu- 
idity always nucleates flrst on the dislocation. A system- 
atic, weakly nonlinear analysis is developed and used to 
derive the one-dimensional Landau theory for superflu- 
idity along the dislocation axis. Finally, we incorporate 
thermal fluctuations into our description of the superflu- 
idity, and determine the effective coupling between two 



"sites" along a single dislocation. The last result moti- 
vates an effective description of the network superfluid- 
ity in terms of a random-bond XY model, and we argue 
that such a model has an ordering temperature expo- 
nentially sensitive to the dislocation density. A series of 
appendices further develops the weakly nonlinear analy- 
sis: Appendix A shows how to improve the naive weakly 
nonlinear analysis using the method of strained coordi- 
nates; Appendix B demonstrates the efficacy of the anal- 
ysis with simple variant of the edge dislocation model; 
and Appendix C extends the equilibrium results to the 
time-dependent Ginzburg-Landau theory. 



II. SUPERFLUID-DISLOCATION COUPLING 
IN LANDAU THEORY 



Following Dorsey, Goldbart, and Toner (DGT),^^ we 
analyze the ordered phase using a Landau theory in 
which the superfluid order parameter V'lx) couples to the 
strain tensor Uij induced by a quenched dislocation. To 
simplify the analysis we assume the solid is isotropic, so 
the superfluid couples to the trace of the strain tensor 
(i.e., divergence of the displacement fleld), which is the 
fractional local volume change of the solid. For a single 
edge dislocation along the z-axis, with Burger's vector b 
along the y axis, the trace of the strain tensor is^^ 



4/i hcosO 
+ X r ' 



(1) 



where /i and A are the Lame elastic constants. We have 
introduced the coordinates x = (r, z), with r in the x ~y 
plane [(r, 6) are polar coordinates in the x — y plane]. The 
Landau free energy functional for the isotropic supersolid 
is then2i'26 



£|V^|2+ l„(r)|^|2+l^|^,|4 



(2) 
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with 



a(r) =ao{to + B 



cosd 



(3) 



Here uq, c, and u are phenomenological parameters (all 
positive); to = {T — To)/To is the reduced temperature, 
with To the mean-field critical temperature for the super- 
solid transition in the absence of the dislocation {to > 
is the normal solid, and to < the supcrsolid); and B 
is a coupling constant (into which we have absorbed the 
elastic constants). Cast in this form, the coupling be- 
tween the dislocation and the superfluid order parameter 
can be thought of as a local change in the critical tem- 
perature due to local changes in the specific volume of 
the solid. 

To simplify the subsequent analysis, we introduce the 
characteristic length scale I = c/aoB, order parameter 
scale X = ao-B/-\/cu, and free energy scale Fq = aoBc/2u, 
and the dimensionless (primed) quantities x' = x/l, 
ip' = ijj/Xi and T = F/Fo'. in terms of the dimension- 
less quantities the free energy becomes 



(4) 



where V{r) = cosO/r, E = —cto/aoB'^, and we have 
dropped the primes on all quantities for clarity of pre- 
sentation. 



III. LINEAR STABILITY ANALYSIS 

Before proceeding further, we analyze the behavior in 
the high-temperature, normal phase. In this phase we 
can neglect the quartic term in the free energy; the re- 
sulting quadratic free energy is 



(fx 



IWI' + [V{r) - E] 



= j (fxii)*{H-E)^, 



(5) 



where the Hermitian linear operator H is given by 

H = -W'^ + V{t). (6) 

We can diagonalize the free energy by introducing a com- 
plete set of orthonormal eigcnfunctions \E'n(x) of iJ, 



E„:^r 



(7) 



where n labels the states, and we assume that the eigen- 
values En are ordered such that Eq < Ei < Equation 

(7) is equivalent to a Schrodinger equation for a par- 
ticle in a potential V{r), and for the dipolar potential 
V{r) — cos9/r the spectrum and eigcnfunctions were 
obtained with extensive numerical work in Ref. 22. Ex- 
panding the order parameter in terms of the eigenfunc- 
tions. 



V'(x) =^A„*„(x), 



(8) 



with expansion coefficients An, and substituting into the 
free energy, after using the orthogonality properties of 
the eigcnfunctions we obtain 



To = Y,{En-E)\An\^. 



(9) 



The free energy is positive as long as En > E, for all 
n. Recall that E = -{c/aoB^){T - To)/To (with ao and 
c both positive), so high temperatures T correspond to 
large, negative values oi E. As we decrease T, E in- 
creases, until eventually we hit a condensation temper- 
ature Tcond at which £'(Tcond) = Eo; below this tem- 
perature the quadratic free energy J"o becomes unstable 
(negative). Rearranging a bit, we have 



Jcond ~ Tq 

To 



-Eo 



aoB^ 



(10) 



If H has negative eigenvalues-i.e., if the equivalent 
Schrodinger equation has bound states-then Tcond > To, 
and the dislocation induces superfluidity above the bulk 
ordering temperature. As emphasized in Refs. 19,21,22, 
the dipolar potential cos6/r, which has an attractive re- 
gion irrespective of the coupling constants, always has a 
negative energy bound state. We are thus lead to the 
surprising and important conclusion^^'^^ that superflu- 
idity first nucleates around the edge dislocation before 
appearing in the bulk of the material. 

Just below the coudeiisution temperature Tcondi the 
nucleated order parameter has the form 



i) = ^o*o(r), 



(11) 



where ^Pq is the normalized ground state wavefunction 
and Ao is an amplitude that is fixed by the nonlinear 
terms in the free energy.^* Substituting Eq. (11) into the 
dimensionless free energy, Eq. (4), we obtain 



T={Eo-E)\Ao\'' + \g\Ao\\ 
where the coupling constant g is given by 

g = j dV*^(r). 



(12) 



(13) 



Minimizing the free energy with respect to Aq, we obtain 

(14) 



Ao = V{E - Eo)/g, 
and the minimum value of the free energy is 

{E - Eof 



^ min — 



25 



(15) 



From the extensive numerical work of Dasbiswas et al?"^ 
we know that for the dipole potential the ground state 
energy is Eo = —0.139 (with the energy of the first ex- 
cited state El = —0.0414), and the coupling constant 
g = 0.0194. 
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To recap-following the work of previous authors, "'^^'^^ 
we have shown that superfluidity always nucleates first 
on edge dislocations, and we have calculated the form of 
the order parameter near i?o, the threshold value of E 
(i.e., at temperatures just below the condensation tem- 
perature). Physically, we imagine a cylindrical tube of 
superfluid, with a radius equal to the transverse correla- 
tion length (of order 1 in our dimensionless units), that 
encircles the dislocation (shown schematically in Fig. 1). 
However, this naive mean-field picture ignores the ther- 
mal fluctuations which destroy the one-dimensional su- 
perfluidity on long length scales. What is needed is an 
effective one-dimensional model for the superfluid, ca- 
pable of capturing nontrivial fluctuation effects. In the 
next section we derive this one-dimensional model using a 
systematic, weakly-nonlinear analysis near the threshold 
So. 




FIG. 1: Schematic showing the dislocation axis (along z) and 
the tubular superfluid region that develops around it. The 
radius of the cylinder is determined by the length scale of the 
ground state wavefunction. The axis of the tube will be offset 
from the dislocation axis. 



IV. WEAKLY NONLINEAR ANALYSIS NEAR 
THRESHOLD 

Within the mean-field Landau theory, the order pa- 
rameter configurations that minimize the free energy are 



solutions to the Euler-Lagrange equation 

— = = -V2^+[y(r)-i?]^+|i^|2V;. (16) 

This nonlinear field equation is difficult to solve, even 
numerically. Instead, we resort to a weakly nonlinear 
analysis^^ near the threshold for the linear instability. 
Our goal is to "integrate out" the modes transverse to 
the dislocation and obtain an effective model for the 
one-dimensional superfluid nucleated along the disloca- 
tion. We then treat thermal fluctuations of this one- 
dimensional superfluid in the next Section. 
We start by introducing a control parameter 

e = E-Eo (17) 

that measures distance from the linear instability. From 
the analysis in the preceding Section, we see that the 
order parameter near threshold scales as e^^^, which sug- 
gests a rescaling of the order parameter 

^ = e'/'0, (18) 

with (f> a quantity whose amplitude is 0(1). Next, note 
that if we had included the plane wave behavior along the 
z-axis in the analysis in Sec. Ill, the coefficient of |A„p in 
the quadratic free energy, Eq. (9), would be -e + P. This 
suggests the important fluctuations along the z-axis oc- 
cur at wavenumbers k ^ e^^^, or at length scales of order 
e~^/^ (i-e., long wavelength fluctuations are important 
close to threshold); the required rescaling is 

z = e-'/\. (19) 

Substituting these variable changes into Eq. (16), and 
writing E = Eq + e, we obtain 

L0 = e[920 + 0-|0|20] , (20) 

where the Hermitian linear operator L is given by 

L = -V'i + V{r)-Eo, (21) 

with the Laplacian in dimensions transverse to z. 
Next, we expand (j) in powers of e, 

= 00 + £01 + e^02 + ■ • ■ ■ (22) 

Collecting terms, we obtain the following hierarchy of 
equations: 



0(1): 


L0O 


= 0, 


Oie): 


L01 


= 9^00 + 00 - 100 P0O 


0{e') : 


L02 


= (1-302)01. 



The solution of the 0(1) equation is the normalized 
ground state eigenfunction, ^o(r); there is an overall in- 
tegration constant Aq{(), so that 

0o = Ao(C)*o(r). (24) 
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Substitute this into the right hand side of the 0{e) equa- 
tion, 



L(j)i = ^od^Ao + *o^o - ^ol^oT^o- 



(25) 



We can determine Aq by left multiplying this cqiiation 
by \l/o, integrating on cPr, and using the fact that L is 
Hermitian, to find 



dlAo + Aa-g\Ao\''Ao = 0, 



(26) 



where g is defined in Eq. (13). This is the solvability 
condition for the existence of nontrivial solutions of the 
0{e) equation. In principle we could solve this equation 
for Aq, substitute back into the right hand side of the 
0{e), and solve the resulting inhomogcncous equation 
to obtain ^2- In practice this is difficult for the dipole 
potential, so we will stop at this level of the perturbation 
theory. 

We can recast Eq. (26) in terms oi z = e~^/^ and 
ip = c^/^Aq as 



aXz) + e<^-.g|(^|V = 0, 



(27) 



which is the Euler-Lagrange equation for the free energy 
functional 



J" 



-I 



dz 



Reinstating the dimensions, we obtain 



-I 



dz 



+ 4I 



(28) 



(29) 



where a = aot, t = {T — Tcond)/7'cond is the reduced 
temperature measured relative to the condensation tem- 
perature, and b = gu. To summarize, we have inte- 
grated out the transverse degrees of freedom (the fluctu- 
ations of which have a nonzero energy gap) and reduced 
the full three-dimensional problem to an effective one- 
dimensional model. In the next section we will study the 
fluctuations of this one-dimensional model, and derive an 
eff'ective phase-only model for a dislocation network. 

A deficiency of our perturbativc approach is the ap- 
pearance of the term (/)o on the right hand side the 0{e) 
equation in Eq. (23). Since (/iq is a zero mode of L, it 
becomes a secular term in subsequent orders of the per- 
turbative expansion. For some systems-e.g., a nonlinear 
oscillator-such secular terms lead to unbounded solutions 
that require regularization using a more sophisticated 



multiple-scale analysis. "'^ In our application the secular 
terms are benign, as the solution decays at infinity (but 
with the wrong decay rate). However, in Appendix A 
we show that the secular terms can be eliminated using 
the method of strained coordinates,^^ with results that, 
on the whole, are equivalent to the results above (but 
with an improved approximation to the order parame- 
ter). In Appendix B we apply our perturbative methods 
to the two-dimensional Coulomb potential ^(r) = — 
for which the calculations can be explicitly worked out 
through O(e^), and we show that the results are in close 
agreement with detailed numerical solutions of the non- 
linear field equation. Finally, in Appendix C we show 
that the same methods may be used to treat dynamic 
equations of motion for the order parameter. 



V. DERIVATION OF THE DISLOCATION 
NETWORK MODEL 

So far we have considered a single edge dislocation, re- 
ducing the full three-dimensional Landau theory to an 
effective one-dimensional theory for a superfluid tube lo- 
calized near the dislocation core. However, a ''He crystal 
will consist of a tangle of dislocations, many of which 
will cross when they come within a transverse correla- 
tion length of each other. Conceptually, we can model 
this as a random lattice (or network) of dislocations, with 
the crossing points serving as lattice sites. While ther- 
mal fiuctuations destroy any long range order in a single, 
one-dimensional tube, the lattice of tubes will generally 
order at a temperature characteristic of the phase stiff- 
ness between adjacent lattice sites. This is the motiva- 
tion behind the models developed by Shevchenko^'' and 
Toner^^. In this Section we revisit the Shevchenko and 
Toner models with a systematic approach to calculating 
the coupling between adjacent lattice sites in the network 
model, and obtain new results on the length scaling of the 
coupling constant. We conclude with some observations 
regarding implications of our results for experiments on 
the putative supersolid phase of *He. 

We start by considering the correlations between two 
points along a single superfluid tube, separated by a dis- 
tance L. Using the notation of the previous section, 
(/?(0) = '-pi and p{L) = (^2 are the values of the superfluid 
order parameter at two sites along the tube. The correla- 
tions are captured by the propagator K{ip2, <^i; L), which 
can be obtained from the functional integral 



K{ip2,ipi\L) = exp{-l3Hes) 




where a, b and c are the parameters of the one- dimensional Landau theory derived in the previous sec- 
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tion, and Hctf is the effective Hamiltonian that charac- 
terizes the couphng between the lattice sites at z = 
and z = L. Cast in this form, we see that the func- 
tional integral for the classical one-dimensional system 
is equivalent to the Feynman path integral for a quan- 
tum particle in a two-dimensional quartic potential (two- 
dimensional because the order parameter (p is complex), 
with z in the classical problem replaced by the imagi- 
nary time T for the quantum system. '"'^ Indeed, previ- 
ous authors have used this analogy to study the effect of 
thermal fluctuations on the resistive transition in one- 
dimensional superconductors, ■^^'■^■^ obtaining essentially 
exact results for the correlation length and thermody- 
namic properties. Consistent with the Mermin- Wagner 
theorem, these authors find no singularities in the ther- 
modynamic properties, and a one-dimensional correla- 
tion length that grows as the temperature is reduced, 
but never diverges. "^'^ These important features are conve- 
niently captured using a simple Hartree approximation,^^ 
in which the qiiartic term is absorbed into the quadratic 
term with the quadratic coeflScient redefined as 



(31) 



where a = ao{T — Tcond)/2cond and ([i^P) a statistical 
average with respect to the effective quadratic theory. 



Carrying out the averaging, we obtain^^ 



(32) 



where ^ — [cjaf-^'^ is the correlation length for the 
one-dimensional system along the superfluid tube, and 
is the cross-sectional dimension of the ID superfluid 
region (the transverse correlation length, as shown in 
Fig. 1). Inserting this result into the Hartree expression 
of Eq. (31), we obtain a cubic equation for the correlation 
length. 



1 



(33) 



where = (c/a)"'^^ is the Gaussian correlation length. 
While ^0 diverges at Tcond) C remains finite at all temper- 
atures (growing as the temperature is lowered), reflect- 
ing the lack of long-range order in the one-dimensional 
superfluid tube. 

Continuing with the Hartree approximation, we can 
find the explicit form of the propagator by exploiting an 
analogy with the partition function for a two-dimensional 
quantum harmonic oscillator; the result is^^ 



ii'(<P2, ;-£') = A:(L)exp < - 



2^sinh(i/0 



b2|' + blP) COSh(X/0 - 2|^2|bl| C0S(^1 - ^2) 
I 



(34) 



where the prefactor k(L) = /3c/2n^sinh{L/^), and ipi = 
|<^i|e*^i and (p2 = |<P2|e'^^. The effective Hamiltonian (up 
to an additive constant) is given by 

H,fi = ^ coth(i/0(l¥'2|' + kiH - Ji2{L) cos{e, - 92), 

(35) 



where 



Jl2(i) 



c\V2\Wl\ 

^sinh(L/^) 
c\ip2\Wi\lL, L/^<1; 
(2c|^2||<^i|/0e-^/^ 



(36) 



The last term in Eq. (35) is the one of interest, as it cou- 
ples the phases at the neighboring sites through an ef- 
fective "ferromagnetic" coupling constant J12 > 0. The 
behavior of J12 as a function of the inter-site separation L 
is one of our important results for small i/^, J12 ^ 1/i, 
reproducing the result of Toner, whereas for large L/^, 
J12 ~ e~^l^. Since ^ is finite for all T, for a sufficiently 
dilute network of dislocations we will always satisfy the 
latter condition-i.e., dilute networks of dislocations pos- 



sess coupling constants exponentially small in the dislo- 
cation density. 

We emphasize that a one-dimensional system with a 

contimioiis symmetry does not exhibit long range order- 
order parameter correlations decay exponentially on the 
scale of the correlation length ^. As a result, we can- 
not replace the full Landau functional with a phase-only 
model (in which one assumes a well-formed order param- 
eter amplitude) , and we need to treat the phase and am- 
plitude fluctuations on the same footing, in the spirit of 
the work by Scalapino et al?^ . In this respect our results 
differ from Toner's^"'^, who used a phase-only treatment 
and found a coupling constant that scales as 1/L. Toner's 
result does apply at short length scales (L <C ^), where 
there is local supcrfliiid order and a phase-only approxi- 
mation can be used. On the other hand, at long length 
scales (L » ^) the exponential decay of correlations re- 
sults in a coupling constant that is exponentially small 
in L. The result in Eq. (36)correctly captures both the 
small and large distance behavior of the coupling con- 
stant. 

So far we have systematically derived an effective 
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Hamiltonian that describes the phase couphng between 
two points (lattice sites) on a single dislocation. To make 
the conceptual leap to the dislocation network, we pro- 
pose that the appropriate model for the network is a ran- 
dom bond XY model of the form 

H=-J2JijCos{ei-ej), (37) 

where (ij) represents nearest-neighbor lattice sites, and 
Jij is a (positive) coupling between the sites that scales 
as e~^iil^ for a sufficiently dilute network of dislocations. 
As noted by Toner, the randomness is irrelevant in the 
renormalization-group sense, ^® and in three dimensions 
we expect the superfluidity in the network to order when 
the temperature is of order the typical coupling strength 
[Jjj]; i.e., fcsTc = 0([Jy]). Again, for a dilute network 
of dislocations (with areal dislocation density XjL?), 
we would find T,, ^ e"-'^/^" ^ a remarkably sensitive 
dependence on the dislocation density. 

How do these theoretical findings compare with the 
experiments? The emerging consensus is that the puta- 
tive supersolid response depends on sample quality and 
preparation (for instance, see Rittner and Reppy'^). As 
an example, the NCRI fraction shows a sensitive de- 
pendence on the sample quality (varying from as much 
as 20% for samples prepared by the blocked capillary 
method and thus having more disorder to 0.5% for those 
prepared under constant pressure,^), while the onset tem- 
perature itself shows a weak dependence on disorder. 
This could suggest that existing experiments arc in the 
high dislocation density regime (L <C 5) where the cou- 
pling constant (and therefore the critical temperature) 
scale as Unfortunately, there have been no sys- 

tematic studies that correlate the dislocation density of 
solid ^Hc samples with the onset temperature for super- 
solid signal. Ultrasound attenuation experiments on solid 
4jjg37,38 suggested dislocation densities ranging from 10''- 
10^ cm^^; at present the dislocation density is not even 
known to within an order of magnitude, which makes it 
difficult to predict from our theory. Moreover, the 
Landau theory parameters for solid "^Hc arc unknown 
(unlike the situation in conventional superconductors, for 
example), compounding the difficulties in directly com- 
paring our results to the existing experiments. However, 
the dislocation network model still serves as a useful con- 
ceptual starting point for understanding some aspects 
of dislocations on Bose crystals, and our work is im- 
portant in establishing the equivalence between the net- 
work model and a more fundamental Landau theory. Fi- 
nally, we note that our theory does not include dynamics 
(however, see Appendix C for a rudimentary discussion 
of dynamics) and we cannot comment on experiments 
that study the rate of superfiow, such as by Ray and 
Hallock^^. However we believe that the model proposed 
above serves as a useful starting point to numerically 
simulate and understand some of these effects. 



VI. SUMMARY 

We have constructed a model for a supersolid based 
on superfluidity induced along a network of disloca- 
tions. Starting from a Ginzburg-Landau theory for the 
bulk solid, we are able to systematically derive a one- 
dimensional equation describing superfluidity along a sin- 
gle dislocation. Then we complete the picture by consid- 
ering a network of dislocations and the effect of over- 
lap of these strands of ID superfluid, that give us back 
superfluid behavior for the bulk. In passing, we note 
the same effect should be observable in superconductors, 
where it should be possible to calculate the various con- 
stants in the Landau theory in terms of microscopic quan- 
tities (something we cannot do for solid ^Hc) . One of our 
more striking results is the sensitive dependence of the 
transition temperature for the dislocation network on the 
dislocation density (i.e., the sample quality). 
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Appendix A: Method of strained coordinates 

In this Appendix we show how to improve upon the 
naive perturbation theory developed in Sec. IV, using 
the method of strained coordinates. We start with the 
Euler-Lagrange equation [see Eq. (16)] 

-\l'^(t) + V{yi)(f)- E(f) + e\(f)\'^(j) = {), (Al) 

where E = Eq+c, with E^ the lowest eigenvalue of — V^-t- 
y(x) (the ground state energy). Assuming V ^ as 
|x| — oo, the bound solutions of this equation behave as 
(j) ~ exp(— V— -Elxj) for |x| — oo; there is a characteristic 
length scale l/^/—E. To account for this scale, define 
new coordinates Xi ~ Xi/l, where I = 1/ ^/T+T/Eq. In 
terms of these new coordinates, we have 

- (1 + e/Eo)Vjc<P + V{IX)(P -E(t> + e\<l)\^(t> = 0. (A2) 

For small e the length scale is / = 1 — e/2£'o -|- 
(3/8)(e/£;o)^ -I- . . .; the potential is then 

y(/X)«F(X)--^X-VxV + 0(e2). (A3) 

Introducing the linear operator 

L = -W\ + V^(X) - Eq, (A4) 
and rearranging things a bit, we have 

U = ^ (^X -VV + V^ 0_ e|</.|20 + O(e2). 

(A5) 
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Now expand cp: 

Substituting into Eq. (A5), we obtain: 



(A6) 



0{1) : Uo 
0{e) : Lcj)i 



0, 



(A7) 



1 /I 



Notice that the right hand sides of both terms do not 
contain by itsclf-there are no secular terms, unhke 
the ''naive" version of the perturbation theory. 

The solution to the 0(1) equation is (pQ = A^q^K), 
where *o is the normalized ground state eigenfunction 
and A is an amplitude. Substituting this into the 0{e) 
equation, 

i</'i = A-^ Qx-VF + y^ *o->l'*o- (A8) 

Now left multiply by ^'o, and use the fact that L is Her- 
mitian: 

(*o, U{) = = A^^ (^i(X • VV)o + {V)o^ - gA^ 

(A9) 



where 



(AlO) 



Note that X ■ W = -X • F, with F the force as- 
sociated with the potential V. The quantum version 

of the virial theorem^^ states that (X • F)o = -2(T)o, 
with T the kinetic energy; therefore, the quantity inside 
the parentheses in Eq. (A9) is the ground state energy 
£^0 = {T)o + {V)o, and we have simply 



= A-gA\ 

orA= 1/^. 

Pulling all of the results together, we have 



(All) 



(A12) 



where 1 = 1/ 1/I + e/Eo = ^/ Eq /E. The strained coor- 
dinate calculation introduces an e dependence into the 
order parameter itself (the "naive" perturbation theory 
gave the same result, but with I = 1), which guarantees 
the correct asymptotic behavior of the order parameter. 



Appendix B: Analysis of a Landau model with a 1/r 
potential 

In this Appendix we solve a simplified version of the 
dipole potential, replacing cos6/r by the attractive two- 
dimensional Coulomb potential V^(r) = —1/r. These two 
potentials share the same length scaling; however, the 



Coulomb potential is rotationally symmetric and the lin- 
ear problem can be solved exactly. The details of the 
pertm-bation calculation follow the general scheme devel- 
oped in Sec. IV. We compare the results of the pertur- 
bation theory with numerical solutions of the nonlinear 
field equation, and find close agreement for a wide range 
of e. 

The energy eigenvalues for the Hamiltonian H = 
— 1/r are given by^^'^^ 



En = — 



(2n+l)2' 



n = 0,l,2,...^ 



(Bl) 



so Eo = —1, with a ground state eigenfunction "^oir) = 
^y2/^^e~'~. The related linear operator [see Eq. (21)] is 



-vi-- + i. 

r 



Expanding </> as before. 



+ 



we have 



Lcpo = 0, 



with the solution 



io(r) =Ao(- 



1/2 



(B2) 



(B3) 



(B4) 



(B5) 



(note that we will ignore the z-dependence in this Ap- 
pendix). Substituting into the left hand side of the 0{e) 
equation. 



Lcpi = (j)Q- |(^o| 00 



3/2 



e-^-". (B6) 



We left multiply this equation by ''' and integrate 

on d^r to obtain = Aq - (l/27r)Ai], so that Aq = \/2tt. 
Substituting back into Eq.( B6), and assuming (pi has 
cylindrical symmetry, we obtain an inhomogeneous equa- 
tion for (pi: 



ld_ 

r dr 



dr 



-4>i + 4>i 
r 



2e-'' - 8e 



-3r 



(B7) 



The explicit solution of this equation (that decays to 
for large r) is 



(t)i{r) = ce-''+e-'^''+re-''+-e 



/•oo t 

ln(2r)+/ — 



dt 



(B8) 

where c is an integration constant. We substitute 0i(r) 
into the right hand side of the O(e^) equation to obtain 



= {i-i(l>l)(pi. 



(B9) 
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We again left multiply by \f2pKe integrate on d^r and 
use the fact that L is Hermitian to obtain the solvability 
condition for the constant c, with the result 



-§ + i+ln(2)-^ln(3) 



'0.75887, (BIO) 



where 7 — 0.577210 is the Euler-Mascheroni constant. 
We have now explicitly calculated two terms in the per- 
turbation expansion; both terms are nonzero at the ori- 
gin, with 00 (0) = 2 and 0i(O) = + ln(2) - | ln(3) = 
-0.047478. 

We have solved the nonlinear field equation for the 
— l/r potential numerically for a wide range of values of 
e = E — Eq = E + 1, using the shooting method.'*'^ The 
results for the order parameter are presented in Fig. 2 
for two different values of e. The two-term perturbation 
theory gives an excellent approximation even for a fairly 
large value of e = 0.43. 




shooting method 
2 term expansion 
1 term expansion 



e=0.09 



10 



shooting method 
2 term expansion 
1 term expansion 



e=0.43 




the origin: 



ei/20o(o)-f £3/201 (0) 
2e^/2 _ 0.0474786^/2 -| 



(Bll) 







T 1 1 1 1 1 1 r 



This result is plotted in Fig. 3, along with our numerical 
results. Again, we see the excellent agreement between 
the two-term perturbation theory and the numerical re- 
sults, even for relatively large values of e. 



2 
1.8 
1.6 
1.4 
1.2 

1 

0.8 
0.6 
0.4 

0.2 



J L 



shooting method + 

2 term expansion 

1 term expansion 



I I I ■ I I I 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

e 



1 IHifjjrT 1 1 1 1 1 1 r 



0.995 

0.99 
0.985 

0.98 
0.975 

0.97 
0.965 

0.96 



+ 



shooting method + 
2 term expansion 



I I I I I I I 



J + 
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FIG. 3: (Color online) Order parameter amplitude ^(0) as a 
function of e. The crosses are the numerical results, the green 
line is the one-term perturbative result, and the red line is 
the two-term perturbative result. The upper plot compares 
the numerical results with the perturbation theory for a wide 
range of e; the two-term perturbative result provides an ex- 
cellent approximation even for values of e as large as 0.8. The 
lower panel is a plot of i/'(0)/2e^/^ as a function of e, which 
highlights the role of the second order term in the expansion. 



FIG. 2: (Color online) Order parameter ^(r) versus r. The 
blue line is the numerical solution, the green line is the one- 
term perturbative result, and the red line is the two-term 
perturbative result. The upper panel compares the results 
for e — 0.09, and the lower panel for e = 0.43. The two 
term expansion provides an excellent approximation to the 
numerical result, even for e — 0.43. 



To get a sense of the efficacy of the perturbation theory, 
we can calculate the amplitude of the order parameter at 



While the naive perturbation theory does an excellent 
job in capturing the overall amplitude of the order pa- 
rameter, it produces the wrong asymptotic behavior of 
the order parameter, as discussed in Appendix A. This 
deficiency is remedied using the method of strained co- 
ordinates, and for the sake of completeness we reana- 
lyze the attractive Coulomb potential problem follow- 
ing the procedure outlined in Appendix A. As noted 
previously, Eq = —1, so the characteristic length scale 
Z = (1 - e)-i/2 = 1 - (l/2)e + (3/8)e2 + . . .. Changing 
coordinates to Xi — Xi/l, and collecting terms, we have 
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at 0(1) 



where 



L(t>o = 0, 



1 
R 



1 = - 



1 d 



R- 



RdR\ dR 



R 



(B12) 



1. (B13) 



The solution is 0o(-R) = At 0(e), we have 

1 



i0i 



2R 



(B14) 



Substituting 0o into the right hand side, left multiplying 
by ^f2pKe~^' ^ and integrating on d^i?, we obtain Aq = 
\/27r, so 



The Oit) equation is then 



-—\R 

RdR\ dR 



R' 



R 



8e 



-3B. 



(B15) 



(B16) 



The explicit solution that decays to zero as — ^ oo is 
[one can show that integration constant c is given by 
Eq. (BIO) as before] 



1 



ln(2i?)+ / dt 



MR) - . - . . . , , 

(B17) 

This is almost the same result as for the naive pertur- 
bation theory obtained in Eq. (B8), where R = \fT^^tr. 
The difference is the third term re"'' in Eq. (B8); this 
is the secular term that is generated in the naive per- 
turbation theory. The strained coordinates method has 
subsumed this term into the first order term; i.e., 



2e" 



2e 



2e 



ere 



(B18) 



To get a sense of the improvement in capturing the cor- 
rect asymptotic behavior of the order parameter, in Fig. 4 
we plot the numerical results obtained using the shoot- 
ing method against the naive and strained coordinate 
perturbation theory. The strained coordinate result is 
indistinguishable from the numerical result in this log 
plot; the naive perturbation theory result clearly decays 
too rapidly for large r. 



Appendix C: Analysis of a time-dependent model 

In this Appendix we generalize the results of Sec. IV to 
derive a one-dimensional dynamical model for the super- 
fluid. For simplicity, we'll assume that there are no con- 
served densities, so the dynamics are described by model 
A (often referred to as time-dependent Ginzburg-Landau 
theory) in the Hohenberg-Halperin classification."''*'"'^ 
The relaxational equation of motion is 



'dt 



-ro 



SF 



?7(x,t), 



(CI) 



10^ 



10"' 



10"- 



10"' 



"T 1 1 1 1 1 1 1 r 



£ = 0.43 



shooting method 
naive 2 term expansion 
s.c. 2 term expansion 



10 



FIG. 4: (Color online) Comparison of the numerical shooting 
method (blue), naive perturbation theory (red dotted), and 
strained coordinates perturbation theory (green dashed) for 
the order parameter 'ip{r) for e = 0.43. The log scale highlights 
the asymptotic behavior of the order parameter-the numerical 
and strained coordinates calculations are indistinguishable for 
large r, while the naive perturbative result clearly decays too 
rapidly. 



where Fq is a relaxation rate and t] is the fiuctuat- 
ing noise term with Gaussian white noise correlations, 
i.e., {r]{x,t)) = and {r]{x,t)*r]{x' ,t')) = 2kBTToS{x- 
x')S{t — t'). As before, to facilitate the reduction to a 
one-dimensional model we introduce dimensionless quan- 
tities, with a time scale i = 2c/{ValB'^) = 2P/{Tc). In 
terms of the dimensionless variables, 



'dt 



where the noise correlations are given by 

{rf (x, i)?7(x', t')) = 2{kBT/Fo)6{^ - ^)5{t - t'). (C3) 

As before, we introduce the small parameter e — E — Eq, 
with V = z — e^-'/^(^, t — e~^T, and fj — e'^/^Tj, to 

obtain 



(C4) 



We expand (f> in powers of e, and at 0(1) we have L(j)o 
0, the solution of which is 



<t)o = ^o(C,^)*o(r). 



(C5) 



Substituting this into the right hand side of the 0{e) 
equation, left multiplying by (/)o, and then integrating on 
d^r, the solvability condition yields 



drAo = d^Ao + Ao- 9\Ao\'^Aq + 



(C6) 



where ^ is the one-dimensional fiuctuating noise term 
(the three-dimensional term with the transverse dimen- 
sions projected out). 



(C7) 



^(C,t)= / d^rM'o(r)77(r,C,T). 
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Using the fact that is normahzed to one, it is straight- 
forward to show that ^ has Gaussian white noise corre- 
lations. Undoing the e scahngs, we obtain our final one- 
dimensional time-dependent Ginzburg-Landau theory, 

9t<^ = 9> + e^-5|<p|V + <e- (C8) 
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